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We study the fermionic Schwinger effect in two dimensional de Sitter spacetime. To do so we 
first present a method to semiclassically compute the number of pairs created per momentum mode 
for general time dependent fields. In addition the constant electric field is studied in depth. In 
this case, solutions for the Dirac equation can be found and the number of pairs can be computed 
using the standard Bogoliubov method. This result is shown to agree with the semiclassical one in 
the appropriate limit. The solutions are also used to compute the expectation value of the induced 
current. Comparing these results to similar studies for bosons we find that while the results agree 
in the semiclassical limit they do not generally. Especially there is no occurrence of a strong current 
for small electric fields. 
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Introduction 


The Schwinger effect, i.e. pair production by a strong electric field, was first studied by Sauter [T]. It has been 
of interest for research over the years and is still today (see e.g. EHH). While these studies are performed in flat 
Minkowski spacetime the development of quantum field theory in curved spacetimes, as a way to combine quantum 
field theory with Einstein’s theory of gravitation [nms], has opened new possibilities. It has been found that the 
gravitational field gives rise to a pair production in a similar way as the electric field [HMD]. This effect is used 
as an explanation for the generation of primordial cosmic inhomogeneities which in turn can explain the large scale 
structure of the universe. Recently also the combination of the effects of the gravitational and electric field have been 
studied in various frameworks pTVl28j . 

As the flat Minkowski spacetime, ZJ-dimensional de Sitter (dS^j) spacetime has constant scalar curvature and is 
maximally symmetric. It is used to describe the early stages of inflation (see for |3D] reviews) as well as the 
late stage of acceleration of the expansion (see [3T]) in cosmology. In addition its study might shed light on the 
understanding of the quantum nature of spacetime [511I33] (see also [TB] page 340ff.). 

The study of the Schwinger effect in dS spacetime is of interest for various reasons. It can be used to study false 
vacuum decay and bubble nucleation muss] as well as to put constraints on magnetogenesis [3^. It is also used 
to test the ER=EPR conjecture using the AdS/CFT correspondence [34]. In addition it can help to understand the 
connection between several methods of regularization better |35| . Very recently a unified thermal picture of Schwinger 
effect in both dS and anti-de Sitter spacetime and Hawking radiation near Reissner-Nordstrom black holes has been 
proposed [HH] . 

The regularized current of scalar QED has recently been computed in dS 2 using Pauli-Villars subtraction in |25| and 
in dS 4 using adiabatic regularization in [35]. While fermionic pair creation in dS 2 has been studied [331 [31 [37] the 
computation of the current has not been performed. Here we will compute the fermionic current in dS 2 and regularize 
it using adiabatic regularization. 

This paper is organized as follows. In section |l] we derive the equations necessary to study the Schwinger effect in 
dS 2 spacetime. We then present a method to compute the pair creation rate semiclassically for general electric fields 
in dS 2 spacetime in section [IT] After this we concentrate on the constant electric field and compute the pair creation 
rate as well as the expectation value of the induced current in section jlll] We conclude in section W To make the 
main ideas clearer we relegated some of the technical calculations to the appendices. 


I. THE DIRAC EQUATION IN dS SPACETIME 

In this section we define the framework in which we will work in this paper. To do so we review the basic equations 
of QED in curved spacetime in EZl Subsequently we will derive the Dirac equation in dS 2 in section jlBj For more 
information on QFT in curved spactime see [131 [T5] . 
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A. QED in curved spacetime 

We start from the action of QED in curved spacetime coupled to a spinor field 
S = [ d^xy/-g{x) 


--R{x) + ^ \^ix)j_^V^ip{x) - (y^-ij^ix)) j^ipix)] - mip(x)'ip{x) - ^Ff,^{x)F^"''{x) 


( 1 ) 


The scalar density of weight 1/2 which generalizes the flat spacetime Lagrangian is 5 := / 6^xC. For the computation 
of the Schwinger effect we assume a background gravitational field defined through the metric as well as a 
background electric field given by 


F^i/{x^ .— d^Ai/{x^ ditAjj^ix'), 

We also introduce the tetrad field e“(a;) with the help of the two dimensional Minkowski metric gab 


gf,^ix) = el{x)el{x)gab- 


The covariant derivative for fermions is given by 


y f,-.= h[d^- -w/ aab ) + ieA^{x), 


where we defined the commutator of the gamma matrices 


<^ab ■■= l^[la,lb] 


as well as the spin connection 


1 


< - 4 [e'’“(x)a^e“ (a;) - e““(a:)a^e^(x) + e““(x)a„e;:(a;) - 

+e'’'"(x)e“^(x)ec^(x)aAe/(x) - e°-''{x)e^^{x)ec^,{,x)dxel{x)] . 

The gamma matrices in curved spacetime 7 ^ are related to the usual ones in Minkowski spacetime via 


7 “ =: 7^e“. 


The canonical momentum is given by 


7r(x) := 


dC 


= V-9ix) ifitp{x)j°. 


( 2 ) 


(3) 


(4) 


(5) 


( 6 ) 


(7) 


( 8 ) 


didoipix)) 

In curved spacetime, the Hermitian adjoint is given by 

'0(x) = il)\x)'-f°. (9) 

The Dirac equation can be derived from the action Q by varying with respect to the field ^/’(x), which gives 

— m) 'ijj{x) = 0 . ( 10 ) 


B. The Dirac equation in dS2 

Two dimensional dS spacetime is described by 

ds^ = — dx^), 


( 11 ) 


^ This definition using the gamma matrix 7 ® of flat spacetime ensures that 'ip{x)ip{x) transforms as a scalar, is real and that the probability 
current = 'ip{x)'y^'il;(x) is conserved. 
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where the scale factor a{r]) depends on the conformal time t] and the constant Hubble factor H following 




(—00 < 77 < 0 ). 


Following from (§ and ( [TT| ) the tetrads for dS space are given by 

= a{r])5l. 

The non-zero components of the spin connection § are found to be 

20 ( 77 )' 


= -wr(x) = 


( 12 ) 


(13) 


(14) 


where prime denotes derivative with respect to tj. 

Using the explicit form of the tetrads (131 and the spin connection (14) the Dirac equation of dS 2 spacetime can be 
derived from ( 10 ) as 


i ( + ^^^7° + ieA^{x)'f ) - m 


ip{x) = 0. 


We choose to work in the Weyl basis 


7° = 


0 1 
1 0 


71 = 


0 r 

-1 0 / ■ 


Using a momentum mode decomposition of the form^ 


® U 2 (r?) 


as well as a solely time dependent electric field A^{x) = (0, 7 l(? 7 )) we find that (15) takes the form 

ih-tpi {?]) + p{T])'ilji{r]) + - ma{r])il^2{v) =0, 

2 0 ( 77 ) 

^ hj q! ij]) 

ih'tfj2{v) -Piv)'>p2{v) + o“rTV'2(??) - ma{rj)^i(rj) =0, 

2 0 ( 77 ) 

where 

^( 77 ) := fc + eA(rj) 

is the kinetical momentum. Decoupling these equations leads to 

>'( 77 ) a'( 77 )' 


ipiiv) + uJkiv) - i^P(7) 


+ ^k{rif + ifipiv) 


piif) a{ri) 

P'iv) _ a'iv) 

p{ri) aij]) 


where we defined the the effective frequency 



[a"(77) 3a'(77)n 


20 ( 77 ) 4a(?7)^ 

+ h^ 

\a"{p) 30077)^ 


20 ( 77 ) 4o(?7)^ 

w?a{'qY. 


■ 01 ( 77 ) =0, 
02 (?7) =0, 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 


The exact form of the momentum decomposition is specified later in J57|l 
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II. SEMICLASSICAL SADDLEPOINT METHOD 

In this section we compute the pair creation rate of general electric helds in dS 2 . This calculation is strongly 
inspired by the one proposed in j39j, where the semiclassical number of pairs is computed in dS 4 . It was remarked 
there that there is an interesting parallel between one-component fields in curved spacetime and two-component helds 
in hat spacetime, for which a generalization of the well known semiclassical techniques (see e.g. gMiiiiiiiaiini) was 
found in [4T]. This comes from the fact that both set-ups have two degrees of freedom. To make this parallel more 
obvious we now consider the held 


4'(x) = ^a{r])tp{x). 


(24) 


This can be seen as an equivalent of the Mukhanov-Sasaki variable in inhation models m- For this held the Dirac 
equation (18)-(191 becomes an equation which is comparable to a Dirac equation in hat space-time with time dependent 


mass 


i{r]) - ma{T])^2iri) = 0 , 
ih'i>2{ri) — piri)'^2{v) ~ rna{T])"^i{ri) = 0 . 


(25) 

(26) 


To obtain the number of created pairs we will now reformulate this equations as equations of the mode functions a{ri), 
Pip). To compute |/?(ry —>■ 0)p we will then perform a multiple integral iteration. Subsequently we will simplify these 
integrals using a a semiclassical saddle point approximation. 

Inspired by similarities between Eqs. (11)-(12) of [H] and (25)-(26), we propose the following ansatz 


= /n ^ / X (a(?7)\/w/c(??) -p(r?)e“ 
^y2u}kip) ^ 



(27) 


(28) 


with the integral 


2 n 

K{p) = - J UJk{T)dT. 


(29) 


The Dirac equation (25)-(26) leads to coupled differential equations for the mode functions by inserting the ansatz 
@-(l2|). They read 


with 



(30) 


(31) 

r(n] - ^k{v)p'{v) 

ma{p) ma{p)uj'i.{p)' 

(32) 


which can be seen as fermionic corrections to the analog bosonic case. 

The initial conditions for the mode functions are chosen such that at past inhnity there are only negative frequency 
modes, i.e. 


«(-oo) = 1, 


P{—oo) = 0. 


It is now possible to iteratively integrate (30) and (31), which leads to 

rO 


Pio) = j2{-i) 




m—0 


2wfe(?7o) 


n 


rVn-l 


Wfc('rn) 




dpr 


^kiVn) 

2Wfc(??„) 


G{pn)e 


— iK{r}n 


(33) 


(34) 
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One can use the fact these integrals are dominated by the classical turning points [ 12 ], which are given by 

‘^kiVp) = 0. (35) 

It is possible to show that for one pair of simple turning points the momentum spectrum of the pair creation rate in 
a semiclassical saddlepoint approximation is given by (see [UHlS]) 


Uk = lim |/ 3 (? 7 )|^ 

77 —>-0 




(36) 


The detailed intermediate steps of the derivation can be found in Eqs. (32)-(38) of |1T|. Observe that in the 

integration contour is closed in the upper imaginary half plane whereas, we (as in [39] 1 close it in the lower imaginary 
half plane because of opposite convention for the phases in (34). 

The above result represents the semiclassical number of pairs created in each mode k for general electric fields in dS 2 . 
In the following sections we will concentrate on the constant electric field, where more explicit computations can be 
performed. 


III. CONSTANT ELECTRICAL FIELD IN dSa SPACETIME 


In this section we will study the Schwinger effect and the expectation value of the induced current in the vacuum 
at negative infinity. We first use the general method of the previous section to compute the semiclassical number of 
produced pairs per momentum mode k for this field configuration in section [III A[ We then construct positive and 
negative frequency mode solutions of the Dirac equation at past and future infinity in section jlll B These solutions are 
used to compute the number of produced pairs per momentum mode k using the Bogoliubov method in section [ill C| 
We then compute the pair creation rate from the previous results in section mE before computing and regularizing 
the current in llll El 

The constant electric field in dS space is described by 


A{r)) = 


E 


since a comoving observer with a four velocity of would measure the field 


= u‘'E^„ = aES^^, 


which has constant field strength Ef^E^- = E'^. 
We introduce 

A = ^ 


7 = 


m 

H 


(37) 

(38) 

(39) 


to characterize the electrical an gravitational field strength respectively. The effective frequency (23) takes the form 

oJkivf = +‘^■k-+ ^ . (40) 

77 77 ^ 

Since we will later construct the positive and negative frequency solutions at past and future infinity we will here 
shortly comment on which form they should take. There form is given by the WKB solution 






UJk{T)dl 


(41) 


Using the explicit form of the effective frequency given in (40) we find that the asymptotic behavior of the positive 


(42) 


and negative frequency modes at past infinity, i.e. for 77 —> —00 is given by 

In the same way the behavior of positive and negative frequency modes at future infinity, i.e. at 77 —> 0, is given by 


i’outiv) 


77 -- 


V\iA ’ 


(43) 


where we defined 


:= + A2. 


(44) 
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A. Semiclassical number of created pairs 


Here we use the semiclassical saddlepoint method described in section |TT] to compute the number of pairs created 
by the constant electric held dehned in (|^. We start by computing the turning points (351, which are given by 


± ^ , 17 

r]p = -- ± —, 


k |fc|’ 


(45) 


where the absolute value of k was introduced such that always denotes the turning point in the lower imaginary 
plane. We now hnd that the real part of is only negative, i.e. inside the integration contour which is used for the 


approximation of (34), if A: and A have the same sign. In |25j this is called pair production in “screening” direction 
because the created pairs would reduce the electric held if we would allow backreaction. In the language of false 
vacuum decay this would be “downward tunneling”. We hnd that as in hat spacetime in the semiclassical limit no 
pairs are produced in “anti-screening” direction which would be connected to “upward” tunneling. This however does 
no longer hold true when we compute the number of pairs using the Bogoliubov method in section [III C| 

Only the imaginary part of K{r]~) is contributing to (36). It is given by 


lm[K{r,;)] = - J (x/ 7 ^ + - |A|) 0(fcA), 


(46) 


where d{x) is the Heaviside step function, 
semiclassical limit is 


Thus we hnd that number of pairs created per mode k (36) in the 


rifc = exp 


-ir(O^-iAi) 


9{kX). 


(47) 


At this point it is insightful to discuss what semiclassical means in this context. The number of created pairs is related 
to the effective action via rik = exp{—S/h). The semiclassical limit is the limit of large action, i.e. \^\ — |A| 3 > 1. A 
necessary condition for this to happen is |/r| 3> 1. Two regimes can be discussed: 

Weak electric field: For A ^ 7 , one hnds S ^ 271(7 ~ I^D- The pairs are mainly produced by gravitational 
(cosmological) pair creation. The hrst term is the usual Boltzmann factor for non-relativistic massive particles at 
the Gibbons-Hawking temperature, while the second term is the correction of the small electric held. The number of 
produced particles gets suppressed by the electric held in this limit. 

Strong electric field: For a strong electric held one hnds S = irmf /(eE), which is the usual semiclassical action for 
the Schwinger effect in hat spacetime. The effect of curvature is negligible in this limit. 

However for the constant held we are not limited to semiclassical calculations but can hnd the solutions of the Dirac 
equation explicitly. We will use this in the following to compute the number of pairs with the Bogoliubov method 
and compare it to the semiclassical result. 


B. Solutions for the constant field 


Here we construct the positive and negative frequency solutions at past and future infinity for the constant field. 
Introducing the new variables 






in (21) and (21) for the constant field (37) we find 


(48) 

(49) 

(50) 

Solutions of these equations are the Whittaker functions IF^±i,^(^), Af^j_i_^(z) [44]. We can make the following 
ansatz for two independent solutions of the decoupled Dirac equation 




1 2 
-I — k- 

1 

-k - 

1' 

-l) 

4 

1 2 

T “ 

z 

1 

-k - 

2 

1' 

K + - 

- 

4 

z 

2 

4 ; 


riz) = 


CiW,_.Jz) 


Z)\ 


(51) 
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Since we solved the two decoupled equations separately we lost the information about the coupling. It can be recovered 
by using the solutions (51) in one of the coupled equations (18)-(19). Using the identity (A 8 ) we find 


Cl/, = = -^7^2/3. 


(52) 


We can show with the help of (A3) that with this values for the constants C 2/3 the solutions (51) are orthogonal 

= 0. (53) 


The value of the remaining constants will be found by asking for normalization after quantizing. 

We can now use the limit of the function 1 Uk,^(z) for \z\ —>■ 00 given in to find the behavior of the solutions ( [slj ) 
at past infinity 


lim ^/>“(z) = lim ^exp , 

V^-oo ri^-oo yC2 Z ^ ^ J \ n J 


lim 4’^{z) = lim 


'q—Y — oo 


C 3 (-z)-'=+ 

^_.,_i)exp 


r,^-oo \Ca{-z)-'^-^ 


kr] . 


(54) 

(55) 


Comparing this to the desired asymptotic behavior of the mode functions given in (42) we find that the positive 
frequency solution ip^iv) is given by for k > 0 and by ip^iz) for k < 0. Thus the positive and negative 

frequency solutions in the asymptotic past can be constructed as 


&) = 


/c > 0 
k<0 


^i„( 2 ) = 


’4^^{z) k > 0 
^/>“(z) k <0' 


(56) 


Following from (53) also these solutions are orthogonal. 

We can now construct the spinor field operator by specifying the momentum decomposition (|T 

/ Hk 

— [b{k)ilj+{'n) + d\-k)i;-{r])] . 

We impose canonical anti-commutation relations 

{'ipaixi,T]),Tr{x[,r])} = ih6{xi - x[)Saf3. 

Using the conjugate momentum ([^ in dS spacetime we find that this is equivalent to 

|■^/'a(a:l,^ 7 ),^/>^(0:1,77) I = - x'i)Sa/ 3 - 

This holds true if the creation and annihilation operators follow the anti-commutation relations 

{b{k), 5(fc')^} = {d{k), c?(fc')^} = 2'xh5{k — k'). 
and the mode functions fulfill the Wronskian condition 

V:+(77)7/:+(77)^' -fT/:" (77)7/:" (77)’^ = ^ 1 - 

0 ( 77 ) 


This is true for the solutions (56) if (see App. E 

C2 = C^ = 


I hH 

m 


e 2 


[k sgn.(fc) 


(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


Accordingly we find that the positive and negative frequency solutions at asymptotic past infinity are given by 




IhH 

m 


, W.+iA^iv) J 

f lU_,+ i,_^(2m) ^ 
Uw_,_^^_^(2iv) 


i’iniv) = 


IhH 

W\ 


W_,+ i,_^(- 2 m) 
/5^1U,_.,^(-2m)\ 

I W,^.J- 2 w) i 


A: > 0 


A: < 0 


(63) 























where we introduced v := \k\ri/h. In an analogous way we can construct the positive and negative frequency solutions 
at ?7 —)• 0 as 


1 hH 




where we introduced 







'Poutiv) = 


r := sgn(fc). 


1 HH 



e^' 




Using the limit of the function Mk,^^{z) for z —>■ 0 given by (A2) we find 

/ 2 \M+| 


lim^+,(^) = limy^e^t 


ri^O r,^0 2 y \k\ 





-/i +2 


(64) 


(65) 


( 66 ) 


(67) 


Comparing to (43) one finds that the modes (64) have the right asymptotic behavior. They can also be shown to 
follow the Wronskian condition (61) by performing steps analogous to the ones found in appendix [B 
Observe that the positive and negative frequency also correspond to the particle and antiparticle solution. This can 
be seen by defining the charge conjugate spinor representing the antiparticle of y’(x) as 'ip‘^{x) := ia 2 fp*{x)^ with the 
Pauli matrix 


CT 2 = 


0 -i 
i 0 


( 68 ) 


This is a physical illustration of Feynman’s picture that antiparticles are traveling backwards in time 


C. Number of pairs using the Bogoliubov method 


We now will use the method of Bogoliubov coefficients to compute the number of created pairs per mode k in dS 
spacetime. Similar methods have been used to compute the pair creation rate in time-dependent fields in flat spacetime 
for general D-dimensional fields in [15] , without an electric field for bosons in dS in [17] and for the constant field in 
dS 2 spacetime for fermions and bosons respectively in |21| and |2S] . In |15j the connection of this technique to kinetic 
theory was shown in the bosonic case. 

We are considering cases where there is a vacuum state for the produced particles in the asymptotic future. This 
requires the rate of change of the background to be small in the asymptotic future, that is 


4M 

4iv) 


r/—>-0 


1^1 ^ and 



(69) 


being small. We see that this is the case when |p,| ^ 1. 

To use the the method of Bogoliubov coefficients we use the fact that the positive frequency mode at past infinity is 
connected to the modes at r] ^ 0 through 


= Oik i’tutiv) + Pk ^out(»?)> 
where the Bogoliubov coefficients are normalized as 

IctfcP + l/3fcP = 1- 


(70) 


(71) 


The coefficients can now be found by putting the explicit form of the solutions (|63| an d (|64 |) in (|70|) and using the 
connection between the Whittaker functions IU„,^(z) and given in 


O-k = 






— rn) y/^ + rn 


e 


e* 


i(r-l) 


Pk = 


r( 2 /r) 


. Using (A4) this leads to 

y/2p 


r(/i — rn) y/^ — rn 


e 2 


i{k-+r 


>e 4 


i(r+l) 


(72) 




























9 


We thus find the number of created pairs per mode to be 

nfe = |/3fc|2 


sinh(7r(|^| + ir^)) 
sinh(27r|/i|) 


(73) 


This result can be shown to be equivalent to Eq. (19) of [2l], where it has been derived in an equivalent way. Comparing 
it to the bosonic result in dS 2 of [35] (see Eq. (2.18)) we find that the only difference is a sinh(7r(|/r| — irn)) instead 
of a cosh(7r(|/x| — irw)). However this difference vanishes in the relevant limit \fi\ 3> 1. 


Comparing to the semiclassical result given in (47) the most striking difference is that according to (73) particles 
can also be produced in “anti-screening” direction which corresponds to “upward” tunneling. This was already found 
in the bosonic case |35|. However in the limit, i.e. |/r| ^ 1,A ^ the two results agree, since pair production in 
“anti-screening” direction gets exponentially suppressed. 


D. Pair creation rate from the number of pairs 


One can compute the pair creation rate from the number of pairs per momentum mode k by integrating 


r := 


1 


dk 


V / 27rh 


nk 


(74) 


where V := a{rf)^drj is the unit two volume of the spacetime. Here we use a derivation of the number of pairs created, 
which relies on an estimate for the moment when most of the particles are created. A good estimate for the moment 
when the particles are created is given when the violation of the adiabaticity is maximal, i.e. when the rate of change 
of the frequency Wfc is extremal. One can show that the maximum gives the estimate for the creation time [25] 




\k\ 


(75) 


Using this the fc-integral can be changed into a time integral. Detailed justification of this procedure can be found 
The pair production rate can thus be estimated from (|73[) by 


r 


cosh(^A) 


27r sinh(27r|/j,|) 

We can also compute the physical number density n of produced pairs at the time 77 with the help of 


1 r 

— dra{r)^T = -. 
Kv ) A-00 H 


(76) 


(77) 


As was already found for bosons this is constant, which shows that the particles created via Schwinger and gravitational 
particle creation compensate the effects of the expansion of the spacetime. That means that the fermion population 
is always dominated by the particles created within one Hubble time, in the limit where our approximations hold. 
Performing calculations in analogy with the ones performed in |26j one finds that in the limit of flat space-time, i.e. 
—>■ 0, one recovers the familiar results for Schwinger pair production in Minkowski spacetime (see e.g. [46]) 


1 |eU| / TT 

hm r = --exp - - 

H^o h TT \ ^ |eA| 


(78) 


We will finish this section by computing the vacuum decay rate. This rate was computed from the imaginary part of 
the one-loop effective action by Schwinger [3]. It has been shown that Schwinger’s result agrees with the canonical 
method for the case of pure electric field [i3H5T] as well as in de Sitter space [53]. In our case, the vacuum decay rate 
is defined as 


T = - / — log(l-rife). 
V / 27Th 


(79) 


Expanding the logarithm and changing the fc-integral to a time integral, as above, it is possible to find the following 
expression for the decay rate 




^)Sinh^(^(|/7| + f) 


27rj 

=±i 7=1 


sinh'^(27r|77|) 


( 80 ) 
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Taking H ^ 0 gives the correct flat spacetime expression 



(81) 


E. Computation of the current 


In this section we will compute the expectation value of the current in a locally inertial coordinate system given by 
= 7 ^. This coordinates were introduced to make the probability density positive semi-definite at each spacetime 
point |3811531154| . The expectation value of the current with respect to vacuum at past infinity is then given by 


J^ = -"(0| [V-(x),7"V'(^)] |0) 


/. 


— oo 
oo 


dk 

2TTh 




-| / ^ - IV'i (»7)P + 1^2 iv)\‘^] ■ 


(82) 

(83) 

(84) 


where ip^{r]) and '0^(7) are the first and second component of tjj^ respectively. Using the diagonal elements of the 
Wronskian condition (|^) we find a connection between the absolute square of the positive and negative frequency 
modes which can be used to simplify the current to 


= -e [ ^ {lipt (»7)P - (vW 


(85) 


Using the explicit form of the positive and negative frequency modes (63) this can be computed as (see Appendix [C|) 

/ 

We will now regularize this current using adiabatic subtraction. 


eH, ( sin( 27 rK) 

J = -1 V - « 

TT \ sin( 27 r/r) 


1. Adiabatic Regularization 


Various methods exist to regularize and renormalize physical quantities. To name some of them, there are proper¬ 
time regularization, dimensional regularization, zeta-function regularization, Pauli-Villars subtraction, point splitting 
regularization (in particular by the Hadamar method) and adiabatic regularization (or subtraction). We propose to 
implement the last one for our problem as it has been done in order regularize the current in the bosonic case in dS 4 
[26] as well as for fermions in flat Minkowski space in [55] . 

Adiabatic regularization was first introduced by Parker to cure the UV divergence and the rapid oscillation of the 
particle number operator [55] . Parker and Fulling generalized it to take care of the UV divergences of the stress energy 
tensor of scalar fields in homogeneous cosmological backgrounds [KTHKH] . 

To regularize the current using adiabatic subtraction we compute the current for slow background variations and 
subtract it from our result (86). To quantify what is meant by slow varying background more precisely, we introduce 
a dimensionless slowness parameter T by replacing the scale factor 0 ( 77 ) by a family of functions aT{r]) := a{r]/T). 
Observe that in the limit of infinitely slow backgrounds, T —)■ 00 , the derivatives of 0 ( 77 ) will tend to zero since for 
integer n 


d^a{r]/T) ^ 1 

d77” T"' 


(87) 


We define adiabatic orders as powers of T~^. In our problem, it is equivalent to count time derivatives and adiabatic 
orders in a given expression. One of the reasons that we did not set h to one here is that since conformal time 
derivative are always multiplied by H, it is actually possible to “formally” pose T := Ijh and to work in powers of h. 
For our purpose, we will expand our modes up to second adiabatic order and subtract it from the current to regularize 
it. To do so we start from the WKB-like ansatz 


tp+ ( 77 ) = Ni 



exp 





p{t) 

2n{t) 




( 88 ) 
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Observe that in difference to the bosonic case of [ 50 ] we are not using a pure WKB ansatz. It was found that for 
fermions, a new kind of ansatz has to be proposed so that the imaginary part of the decoupled Dirac equation ( 21 ) 
is canceled. This has been used for fermions under the influence of an electric field in flat spacetime in |55j and in 
curved (Friedman-Lemaitre-Robertson-Walker) spacetimes without electric field [35] (see also [SO])- The ansatz (88) 
is a combination of the two previous ones. By putting the ansatz in the decoupled Dirac equation ( 21 ) we find a 
reparametrization of it in terms of D(?7), namely 


11(77)^ — 01(77)^ = h ^ 


20(77) 


, Pijl) \ _ a'ijl? /3 77(77) _ p{r]f \ a'ij]) 77'(77) 

0(77)/ 0(77)2 \^4 211(77) 4 D( 77 ) 2 y 20(77)11(77) 

77"(77) + fl"(77) ^ 30'(77)2 4 


4 


(p'{t) a'{t) 


1 _ 
■p'{iiy 


( 89 ) 


ri(t)2 \p{t) a{t) J 211(77) 40(77)2 

We can now expand ( 89 ) to find 0 ( 1 ) = uj{t) 4 0 {T)~‘^. The (n -I- l)-th order can be found by iteratively using the 


7 i-th order solution on the right hand side of ( 89 ). 


For the second component of the spinor we use the ansatz ( ] 88 | ) in the coupled Dirac equation ( 18 ) to find 

0(77) -k 77(77) 


7/7 j" ( 7 ?) 

7/)+(77) 


7770 ( 77 ) 


iH 

2ma{ri) 


O'(77) +77'(77) 
0(77) 


Q'(7?) 

0(77) 


1 


( 90 ) 


Solving the off-diagonal element of the Wronskian condition ( 61 ) for (77) and using it in one of the diagonal elements 
we can show that 


+ |7/’^(7?)P = 


1 


77 ( 77 )’ 


( 91 ) 


Using this normalization condition we can now write the current ( 85 ) in terms of the fraction ( 90 ) 


r = - 


dk 


1 - 


1^2 






( 92 ) 




To perform the adiabatic expansion we can use the fact that ^{rjY ”<^(77)2 is of second adiabatic order, which follows 
from ( 89 ). It is possible to write 12(77) = \/a;(77)2 -p [11(77)2 — 07(77)2] and then expand the square root to find 


11(77) =^(77) 4 


11 ( 77)2 - w(77)2 


207 ( 77 ) 


0 {T) 


-4 


Using this we can expand ( | 92 [ ) to second adiabatic order 

07 ( 77 ) - 73 ( 77 ) 


r = 



207 ( 77)3 


[ 11 ( 77 ) 2 - 07 ( 77)2 


772 [777a'(77)73(77) — 773,0(77)73'(77)]^ 
¥ 07 ( 77)6 


+ 0 {T) 


-4 


( 93 ) 


( 94 ) 


where for 11(77)2 —uj{ri)‘^ one can use ( 89 ) and replace 11(77) by 07(77) on the right hand side. We can now compute this 
for the constant electric field and find 


= 


a{v) J-c 


ZTTnujyr}) TT n 


-4 


( 95 ) 


For the bosonic case, the same counter term was found using the Pauli-Villars regularization method in dS2 in [ 25 ] . 
The adiabatic regularization for bosons, in dS2, is very similar to the calculation presented above and the final result 
can be found to be the same. 


Performing the adiabatic subtraction for the current ( 86 ) we find 


eg sinh(f A) 
^ ^ sm(27rfi) ' 


( 96 ) 


In the next subsections, we will discuss the properties of this current in detail. Performing the limit of strong electrical 
and strong gravitational field we underline the effect of the respective contributions to the total Schwinger effect. For 
a plot of the current ( [ 9 ^ as a function of A for different values of 7 see Fig. There we also show a comparison to 
the bosonic current of | 25 | which is discussed in section HIE 3 and plotted as dotted lines. 
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2. Strong and weak electric field limit 


If we look at the limit A —>■ oo for fixed 
asymptotic behavior 


7 , we find that the current is dominated by the electric field and has the 


- _ : 

'^reg 


eH\ 
TT h' 


(97) 


which is independent of the mass. This linear behavior can also be seen in Fig. where all the curves align for large 
A. In the bosonic case in dS 2 this linear behavior is also present |25j whereas in dS 4 a quadratic behavior is found 
which leads to a linear behavior for the conductivity (defined as J/A). This can be used to impose strong constrains 
on magnetogenesis scenarios |26j . For fermions in dS 4 the same might appear. The strong field limit is the same limit 
as the small mass limit "f H for which we find 


TT I h 


H f 2 ty 

— — TTCoth —A 
2A \ H 




(i)' 


(98) 


For 7 = 0 we find exactly the linear response which is analog to the asymptotic behavior (971 since for strong enough 


fields the effect of the mass is negligible. The linear response is the usual flat spacetime response. Indeed the Schwinger 
pair creation rate for massless carriers in flat spacetime is given by 


F = 


eE 

2'Kh' 


(99) 


which leads to the induced current 


J = 


e^Et 

ttH 


( 100 ) 


In an expanding spacetime the current can be computed in a comoving frame but will then be diluted in the physical 
frame. Naively making the substitution t —> l/i7, gives the linear response found in (97). This substitution will be 
used again in the context of the flat spacetime limit in section HIE4 This shows that taking the limit 7 = m/i7 —>■ 0 
is equivalent to take the limit H ^ 0 and m = 0 which is the flat spacetime limit for massless particles. This illustrates 
the fact that a massless fermion is conformally invariant. 

The limit of strong electric fields also allows us to compare to the result of the pair creation rate found in section 
|III D[ An approximation of the current for relativistic particles is given by 


J « 2en = 2e—. 

H 


( 101 ) 


Since the particle picture used to derive the pair creation rate (76 1 only holds in the limit |/r| 3> 1 which is equivalent to 
■^ 7 ^ + A^ ^ h and the assumption of relativistic particles only holds for small masses, the limit in which a comparison 
is possible is the strong field limit. 


Using the pair creation rate (76) to compute the estimated current (101) in the limit |p,| ^ 1, A ^ we find 


7 I 

J « -|/r| exp 

TT 


-27r ( ImI - -y 


( 102 ) 


which agrees with the regularized current (96) in this limit. 


Expanding (96) for small electric helds A <C we find 


ei7 27 A / A 

sinh(f7) 


(103) 


As also visible in Fig. [^for small electric fields the mass begins to have an effect on the current. We find that for 
small electric fields the current gets strictly decreased by increasing mass. This is in contrast to what was found in 
the bosonic case of [25] and we will discuss this in the next section. 
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A 
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FIG. 1: Regularized current for fermions and bosons (dotted) given by (96 1 and (1041 respectively, as a function of A for 
different values of 7 . One sees that for small electric fields there is a difference while all curves have the same asymptotic limit 
(971 for negligible mass. Observe that only the bosonic case shows the “infrared hyper-conductivity”, i.e. a large current for 
small electric field and mass. The curve with the linear response is given by 7 = 0 in the fermionic case while in the bosonic 
case it is found to be at 7 = O.Sfi. 


3. Comparison to the bosonic case 


We can compare the regularized current (96) to the scalar one in dS 2 given by 

eH sinh (^A) 


rx 


TT sin (27rcr) 


(104) 


where tr = y /i 

Observe that this is the same as the fermionic current ( |96[) if we replace cr —> /r. This means that the two currents 
agree in the limit |/i| 3> 1 since there cr « /x (see also Fig. Ill), as we have already shown for the pair creation rate in 
section EnDl However we find a different behavior for electric fields small with respect to the Hubble constant. This 
is due to the fact that n ^ a and is in difference to the flat space case where the pair creation of a constant electric 
field is identical for scalar and spinor QED. 

The most striking result of the difference is, that the bosonic current is enhanced for small electric fields and values 
of the parameter 7 < 1/2. In this effect was given the name “infrared hyper-conductivity”. We do not find this 
behavior for fermions. Mathematically this is due to the fact that a in difference to the parameter /x can become real 
for 7 < 1/2. In the same way the linear response for the bosonic case is found for 7 = 1/2 whereas it is found for 
7 = 0 for fermions. Since in flat space-time the linear response arises for massless particles the fermionic result is not 
as peculiar as the bosonic one. 


The difference between a and /x comes from the last term proportional to in (21)-(22). For the constant electric 
field (37) it evaluates to 1 /( 477 ^) and accounts for the 1/4 term in the Whittaker equations (49l-(50). These terms 


are absent in the Klein Gordon equation and thus an additional factor of 1/4 has to be introduced in the variable a. 
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4- Vanishing gravitational field: flat spacetime limit 


For a comparison with the flat spacetime case the parameters A and 7 deflned in (391 are not convenient since they 
diverge in the flat space limit H ^ 0. We therefore define the electric field e in units of the critical electric field as 

(105) 


= nA. 


Performing the the flat spacetime limit if —)■ 0 of (96) we find 


em 


lim = lim 
H^o ® H^o H 


e _ 
e I 


(106) 


This current diverges, as it is expected for the current of an electric field in Minkowski spacetime, which was turned 
on at past infinity. This is due to missing Hubble dilution. One can however compare it to a current of an electric 
field which was turned on at finite time t, as it was done in the bosonic case [55]. Using the substitution 1/H ^ t 
introduced above, one finds agreement with the current induced by the Schwinger effect in flat Minkowski 


IV. CONCLUSION AND PERSPECTIVES 


Recently many generalizations of the original Schwinger effect in de Sitter spacetime were proposed. While most 
of these studies focus on bosonic pair creation, we treat fermionic pair creation in dS 2 in this article. After reviewing 
the basic equations in section |I| we perform a first derivation of the number of produced pairs based on a semiclassical 
saddle-point approximation. This derivation was already proposed in the more general case of dS 4 in [35] and is 
straightforwardly adapted to dS 2 . Observe that the results in section |n| are independent of the specific form of the 
field and can be used to compute the pair creation rate for general time dependent fields. 

In the rest of the paper we concentrate on the effect of a constant electric field. First we use the techniques of section 
jTTjto compute the semiclassical number of pairs per momentum mode k. We then construct the positive and negative 
solutions in asymptotic past and future to use the method of Bogoliubov coefficients to compute the number of created 
pairs again. Here we find a first difference to the bosonic case of [25j in equation (73). However particle creation is 


only well defined for well defined particle states in the asymptotic future. We show that it is the case when |^| ^ 1. 
In this limit the fermionic and bosonic number of created pairs agree with each other as well as with the semiclassical 
result. 

This leads us for our first major conclusion: bosonic and fermionic pair creation are the same in the semiclassical 
approximation. In flat spacetime, this result is well known for one component fields with one pair of turning points 
but has been shown to be false for non constant fields with more than one component mi- It could be interesting 
to explore more complicated electromagnetic field configurations in curved spacetime and check if the equivalence of 
semiclassical bosonic and fermionic pair creation rates still holds true. 

Computing the pair creation rate from the number of pairs created per momentum mode k, a diverging momentum 
integral appears and to avoid the pair creation rate to diverge, we proposed to transform the momentum integral into 
a conformal time integral via a discussion of the non-adiabaticity. This heuristic step is necessary to find a coherent 
result, but is a sign that the pair creation rate per unit two volume is not the best quantity to describe the Schwinger 
effect in curved spacetime. This is due to the fact that the very definition of pair creation time does not exist. 

The quantity which accounts for Schwinger effect beyond the particle picture is the induced fermionic current. It 
is a more relevant quantity than the usual pair creation rate per unit two volume because it is not plagued by the 
absence of a clear definition of pair creation time. The final result is presented in (96). It has been renormalized by 
a adiabatic subtraction. The correction turns out to be the same as the one found in the bosonic case of |5^ where 
the Pauli-Villars method was used. We find that the fermionic regularized current is very similar to the one found in 
the bosonic case. However the effect of “infrared hyper-conductivity” is not present for fermions. 

It would be interesting to investigate the relations between bosons and fermions in 4 dimensions. However we did not 
find an analytic solution of the coupled differential equation arising in this framework. Performing a study analog to 
the one in this work or the bosonic case of |26j seems still possible numerically and is left for future work. The pair 
creation rate derived in the semiclassical saddle-point approximation has already been calculated and shown to agree 


® See e.g. Eq. (5.22) of 1471 or Eq. (2.14) of 1611 for the bosonic current in four dimensional flat spacetime and Eq. (5.32) of 1481 for the 
bosonic current in two dimensional flat spacetime. 
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with the bosonic one in [35] . 

As the gravitational and electrical field were taken to be external, an interesting extension would be to take their 
variation into account, i.e. use the induced fermionic current in the generalized Maxwell equation and the resulting 
number of pairs in the Einstein equation. This could be used to find specific forms of electric fields or specific 
classes of spacetimes which favor or disfavor pair creation. The fact that Pauli-Villars (in the bosonic case [35] 1 
and adiabatic regularization give rise to the same counter term, shows that pair creation is a possible framework to 
compare regularization techniques in curved spacetime. Other interesting issues not treated in this paper but worthy 
to investigate are cosmological applications of our result. The created pairs might account for the asymmetry of 
matter/anti-matter in our universe with a modification of an Affleck-Dine mechanism or a specific model of preheating. 
They could also give hints on the evolution of an accelerated period of expansion and how matter and gravitation 
interact in such periods. 
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Appendix A: Properties of the Whittaker functions 


For the construction of the mode functions we use the Whittaker functions and These well studied 

functions have a number of properties which can be found in the literature (see e.g. [44j 1. The ones we use throughout 
the paper are summarized in this appendix. 

The functions have the limiting behavior for |z| —> oo 


lim VPk,^(z) = lim e for |arg(z)| < -tt, 

\z\—>-00 |2|—>-00 Z 


(Al) 


as well as for z —>■ 0 


limM,, „(z) = lim z^+^Z^. 

Under conjugation they behave as 

[W^A^)]* = lU«.,^.(z*), [M.Az)T = M«.,^.(z*). 

and they have the properties 

1U„ _^(z) = 1 Uk,^(z) = ±ie=^’"‘Af_«,,^(z). 

The function WU_^(z) can be written as an integral in the Mellin-Barnes form 


^ ^ i-ioo 2i7r r (i + - k;) r (i - /X - k) 


The two functions are connected through 


W.A^) = ■ .p = 


The Wronskian is given by 




— ^ — ’ r(i-i-/x —k) 


(e±‘-z) dW.A^) 


dz 


By using r(n -|- 1) = r(n)n in (A5) one can find 

2k -I- 1 — z 






2{A-k^) dz 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 




(AS) 
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Appendix B: Wronskian condition 


In this appendix we show how to derive the value of parameters (621 of the positive and negative frequency 
at asymptotic past (56) so that the Wronkskian condition (61) holds. We therefore define 

M •= + tlj~ {r])^. 


Now we can use the behavior of the Whittaker function under conjugation given by (A3) and the specific foi 
solutions (56) given in (51) and (52) to find 


V'l ivT = '02 (^) • S , c! 


— ^ for fc > 0 


+ ^ for fc < 0 


iv)* = ■ s cl 


+ ^ for fc > 0 


— ^ for fc < 0 


where 'ijj^ (rj) and '0^ ( 77 ) are the first and second component of 00 respectively. This can be used to find 

Mi2* = M21 =0+ (?7)00 (77) ^ ^ 

O 2 O 3 


Now requiring the Wronskian condition (611, i.e. M 12 = M 21 = 0, we find 

|C2P = |C3|0 


Using (B2) for the diagonal elements of M we find 


C? 


Mil = M 22 * = --^sgn{k) [ 0 j^(? 7)02 (?7) -01 (^)0^(?7)] • 
(-/2 


We can now use the Dirac equation (18) to bring this into the form 

ih a 


Ml 1 = M 22 = — - 


1(77) C2 


sgn(fc) 0+ ( 77 ) 0 ^ ( 77 ) - 01 ( 77)00 ( 77 ) 


Using the solutions (56) we find 


Mil = M 22 * = — 


mayrj) C 2 [077 drj 


using the explicit form of the solutions (51) and (52) this becomes 


Mil = M 22 * = —i 


.ICaP 2k 


airf) Hh 


_„(-z) dWi^_i Jz) 

- dz - wl 


dz 


Using the Wronskian (A7) we find 


Mil = Moo* = _ IC'al ‘^\k\ iTrsgn{k)K 

ivin 1 V 122 \^r,) Hh a{d) Hh 


Now requiring the Wronskian condition (61), i.e. Mu = M 22 = 0 ( 77 ) we find 

2 _ i^Ksgn(/c) 

- 2\kf 


Choosing a physically irrelevant phase this leads to (62) using (B4). 
Observe that from (B2) and (|B4|) we find 


= 102 Wit 


solutions 

(Bl) 
of the 

(B2) 

(B3) 

(B4) 

(B5) 

(B 6 ) 

(B7) 

(B 8 ) 

(B9) 

(BIO) 


l0tWt = 101 Wt- 


(Bll) 
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Appendix C: Computation of the integral for the current 


If we change the variable in (85) to v := \k\r]/H we find 

dv 


r = - 


f ) 

Jq ZtT \ k>0 k>0 k<0 k<0/ 


V Jo 27r 

2 e dv 




fc >0 


J, 

fe<0/ 


(Cl) 

(C2) 


where we used the normalization (91) of the modes. Using the specific form of the solutions (63) this takes the form 

..2 H-„ roo , 2 


JX ^ ^ I ^ 

2n Jo V 


r e 


r=-ia 




(C3) 


Using the Mellin-Barnes form of the Whittaker function (A5) we can write the current as 


ZTT ^->00 , 


dv 

V 


r°° ds_ 

—2i7r 


dt 

2m 


{2iv) ®(— 2 iu) * ^ re 

r=-l,l 


r(l/2 + fi + s)r(l/2 - ^ + s)r(l/2 -m-s) r(l/2 - /r + t)r(l/2 + fi + t)T{l/2 + rK-t) 


(C4) 


r(l + /r — rK)r(l — fi — tk) 

We now perform the integral over v 

7 ^ He 


r(l — fi + rK)r(l + /r + rn) 


H =- 


lim 


(27r)3 ^->00 r(l + /r — K)r(l — p. — K)r(l — ^ + ^)r(l + /r + «) 

p'lOO _ 

I ds re"™e"‘5"r(l/2 + ^ + s)r(l/ 2 -/r + s)r(l/ 2 -rK-s) 
1 

(20-«-‘ 


(C5) 


r=-l,l 


dte'^* r(l /2 -n + <)r(l /2 + H + <)r(l /2 + rK-t) ' 


s +1 


We can close the contour of the integral over t in the t > 0 -plane in order to use the residue theorem. The contributing 
poles thus are t = —s and t = —rn -I- 1/2 -|- n. Due to the ^—>-00 limit the integral is only non-zero for t -|- s < 0. Since 
we could close the integral over s in a similar way in the s > 0 -plane for t ^ —s the only pole which gives non-zero 
contribution is t = — s. Using the residue theorem we thus find 


r = - 


7 ^ \eH 


r°° r(l/2 + H + s)r(l/2 -fi + s)r(l/2 - /r - s)r(l/2 + fi-s) 

/_ioo ^ r{i + fi — K)r{i — fj, — k)t{i — ^ + K)r{i + fj, + k) 

r e“‘’^*'’''‘“''®^r(l/2 — rn — s)r(l/2 -|- tk -|- s) 


1 eH. 


X 

T-=-l,l 
/*ioo 


ds 


sin[(/i -I- k ) tt ] sin[(/i — k ) t :\ 


47 r cos[(/r -|- s) 7 r] cos[(/r — s) 7 r] ^ cos[(rK -|- s) 7 r] ’ 


E 


^—i7r{rK,-\-s) 


(C 6 ) 


(C7) 


To solve this integral we now write trigonometric functions as exponentials and change the integration variable to 
X := exp(i 7 rs). This leads to a integral which can be solved using the standard decomposition theorems for rational 
fractions 


r = - 


7 


eH. 


jj? — TT^ 


dX 


sin[27rK] sin[(^ -|- k)tt\ sin[(/r — K)7r]W^ 

(A2 -p e2i’^'=)(X2 -p e-2i’"«)(A2 -p e^'^^^^){X‘^ + ' 


eH. 


sin(27rK) 
TT \^^sin(27r^) 


=- 1 A* 


(C 8 ) 

(C9) 
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